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Energetics and quantized conductance in jellium-modeled 
nanowires are investigated using the local-density-functional- 
based shell correction method, extending our previous study 
of uniform-in-shape wires [C. Yannouleas and U. Landman, 
J. Phys. Chem. B 101, 5780 (1997)] to wires containing 
a variable-shaped constricted region. The energetics of the 
wire (sodium) as a function of the length of the volume- 
conserving, adiabatically shaped constriction, or equivalently 
its minimum width, leads to formation of self-selecting magic 
wire configurations, i.e., a discrete configurational sequence of 
enhanced stability, originating from quantization of the elec- 
tronic spectrum, namely, formation of transverse subbands 
due to the reduced lateral dimensions of the wire. These sub- 
bands are the analogs of shells in finite-size, zero-dimensional 
fermionic systems, such as metal clusters, atomic nuclei, and 
^He clusters, where magic numbers are known to occur. These 
variations in the energy result in oscillations in the force re- 
quired to elongate the wire and are directly correlated with 
the stepwise variations of the conductance of the nanowire 
in units of 2e^//i. The oscillatory patterns in the energetics 
and forces, and the correlated stepwise variation in the con- 
ductance are shown, numerically and through a semiclassical 
analysis, to be dominated by the quantized spectrum of the 
transverse states at the narrowmost part of the constriction 
in the wire. 



PACS numbers: 73.20.Dx, 73.40.Jn 



I. INTRODUCTION 

Understanding of the physical origins and systemat- 
ics underlying the variations of materials properties with 
size, form of aggregation, and dimensionality are some 
of the main challenges in modern materials research, of 
ever increasing importance in the face of the accelerated 
trend towaipdjrniniaturization of electronic and mechani- 
cal devices. Era 

Interestingly, it has emerged that concepts and 
methodologies developed in the context of isolated gas- 
phase clusters and atomic nuclei are often most useful for 
investigations of finite-size solid-state stjpiptures. In par- 
ticular, it has been shown most recentlyO'Q through first- 
principles molecular dynamics simulations that as metal- 
lic (sodium) nanowires are stretched to just a few atoms 
in diameter, the reduced dimensions, increased surface- 
to-volume ratio, and impoverished atomic environment, 
lead to formation of structures, made of the metal atoms 



in the neck, which can be described in terms of those ob- 
served in sroall gas-phase sodium clusters; hence they 
were termedO^cl as supported cluster-derived structures 
(cds). The above prediction of the occurrence of "magic- 
number" cds's in nanowires, due to characteristics of elec- 
tronic cohesion and atomic bonding in such structures of 
reduced dimensions, are directly correlated with the en- 
ergetics of metal clusters, where magic-number sequences 
of cluster sizes, shapes and structural motifs due to elec- 
tronic and/or geometxijc shell effects, have been long pre- 
dicted and observed .LTlI 

These results lead one directly to conclude that other 
properties of nanowires, derived from their energetics, 
may be described using methodologies developed pre- 
viously in tlic context of clusters. Indeed, in a pre- 
vious letter,E3 we showed that certain aspects of the 
mechanical response (i.e., elongation force) and elec- 
tronic transport (e.g., quantized conductance) in metal- 
lic nanowires can be analyzed using the local-density- 
approximation (LDA) -based shell correction method 
(SCM), developetLj^d applied previously in studies 
of metal clusters.lalid Specifically, we showed that in 
a jellium-modelled, volume-conserving, and uniform in 
shape nanowire, variations of the total energy (partic- 
ularly terms associated with electronic subband correc- 
tions) upon elongation of the wire lead to self- selection 
of a sequence of stable "magic" wire configurations 
(MWC's, specified by a sequence of the wire's radii), with 
the force required to elongate the wire from one config- 
uration to the next exhibiting an oscillatory behavior. 
Moreover, we showed that due to the quantized nature 
of electronic states in such wires, the electronic conduc- 
tance varies in a quantized step-wise manner (in units of 
the conductance quantum go = 2e^//i), correlated with 
the transitions between MWC's and the above-mentioned 
force oscillations. 

In this paper, we expand our LDA-based treatment to 
wires of variable shape, that is allowing for a constricted 
region. From this investigation, we conclude that the 
above self-selection principles and the direct correlations 
between the oscillatory patterns in the energetic stability, 
forces, and stepwise variations of the quantized conduc- 
tance maintain for the variable-shaped wire as well, with 
the finding that underlying these oscillatory patterns and 
correlations are the contributions from the narrowmost 
region of the wire. Furthermore, this finding is analyzed 
and corroborated through a semiclassical analysis. 

Prior to introducing the model studied in this pa- 
per, it is appropriate to briefly describe certain previous 
theoretical and experimental investigations, which form 
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the background and motivation for this study. Atom- 
istic descriptions, based on realistic interatomic interac- 
tions, and/or first-principles modelling and simulations 
played an jessential role in discovering the formation of 
nanowires,ll3 and in predicting and elucidating the micro- 
scopic mechanisms underlying their mechanical, spectral, 
electronic and transpori-properties. 

These predictionsEj^EJ [particularly those pertaining 
to generation of nanowires through separation of the con- 
tact between two materials bodies; size-dependent evo- 
lution of the wire's mechanical response to elongation 
transforming from multiple-slips for wider wires to a suc- 
cession of stress accumulation and fast relief stages lead- 
ing to a sequence of structural instabilities and order- 
disorder transformations localized in the neck region 
when its diameter shrinks to about 15 A; consequent 
oscillations of the elongation force and the calculated 
high value of the resolved yield stress ( ~ 4 GPa for Au 
nanowires; which is over an order of magnitude that of 
the bulk), as well as anticipated-jelectronic quantization 
effects on transport propertiestSEj] have been corrobo- 
rated in a number of experimrpts— iising scanning tup-, 
neling and force micros(pipy.EJ't^"Eil break junctions,EHl 
and pin-plate techniquesli3n3 at ambient environments, 
as well as under ultrahigh vacuum and/or cryogenic con- 
ditions. Particularly pertinent to our current study are 
experimental observations of the oscillatory behavior of 
the elongation forces and the correlations between the 
changes in the conductance and the force oscillations; see 
especially the simultaneous measurements of force and 
conductance in gold nanowires in Ref. |2^, where in addi- 
tion the predicted "ideal" value of the critical yield stress 
has also been measured (see also Ref. 21). 

The LDApj^llium-based model introduced in our pre- 
vious papert2l and extended to generalized wire shapes 
herein, while providing an appropriate solution within 
the model's assumptions (see section II), is devoid by 
construction of atomic crystallographic structure and 
does not address issues pertaining to nanowire forma- 
tion methods, atomistic configurations, and mechanical 
response modes [e.g., plastic deformation mechanisms, 
interplanar slip, ordering and disordering mechanisms 
(see detailed descriptions in Refs. p^|JT^ a nd |lj, and a 
discussion of conductaiune dips in Ref. |18|), defects, me- 
chanical reversibilityj^tj and— roughening of the wires's 
morphology during elongationlij] , nor does it consider the 
effects of the above on the|,electron spectrum, transport 
properties, and dynamics.Q Nevertheless, as shown be- 
low, the model offers a useful framework for linking in- 
vestigations of solid-state structures of reduced dimen- 
sions (e.g., nanowires) with methodologies developed in 
cluster physics, as well as highlighting certain nanowire 
phenomena of mesoscopic origins and their analogies to 
clusters. 

In this context, we note that several other treatments 
related to certain of the issues in this paper, but em- 
ploying |&j8e|electron models, have been pursued most 
recentlyBil'Ea In both of these treatments an infinite con- 



fining potential on the surface of the wire is assumed 
and only the contribution from the kinetic energy of the 
electrons to the total energy is considered, neglecting 
the exchange-correlation and Hartree terms, and elec- 
trostatic interactions due to the positive ionic (jellium) 
background. A comprehensive discussion of the limita- 
tions of such free-electron models in the context of calcu- 
lations of electronic structure and energetics (e.g., surface 
energies) of metal surfaces can be found in Ref. |2^. 

In section II. A., we outline the LDA-based Shell Cor- 
rection Method, describe the jellium model for variable- 
shaped nanowires, and derive expressions for the ener- 
getics of such nanowires (density of states, energy, and 
force). Numerical results pertaining to energetics, force, 
and electronic conductance, calculated as a function of 
elongation for variable-shaped sodium nanowires, are 
given in section II. B., including a discussion on the main 
finding that the contribution from the narrowmost part 
of the constriction underlies the properties of these quan- 
tities and the correlations between them. These correla- 
tions between the energetic and transport properties and 
their dependence on the narrowmost part of the nanowire 
are further analyzed in section III, using a semiclassical 
treatment. We summarize our results in section IV. 



II. DENSITY FUNCTIONAL DESCRIPTION OF 
JELLIUM NANOWIRE 

A. Theory 

1. Shape of Constriction 

Consider a jellium nanowire with circular symmetry 
about the axis of the wire (z axis). The wire may contain 
a constricted region (see Fig. 1), that is a section of length 
L where the cross-sectional radius a{z) varies along the 
axis as 

a{z) = ao + (i?o - ao).f{z) , -1/2 <z<L/2 , (1) 

with f{—z) — f{z) (the z = plane passes through the 
middle of the wire) and f{±L/2) = 1. Rq = a(±L/2) is 
the uniform radius outside the constricted section, and 
flo = a(0). In this paper, we take a parabolic shape 
f{z) — (2z/L)^ for the description of the constricted re- 
gion [a wire of uniform cross section throughout corre- 
sponds to f{z) = 1]. 

We also assume that elongation of the wire occurs in 
the constricted region while maintaining its volume con- 
stant (this is supported by MD simulations), namely by 
requiring that 



L/2 



a^{z)dz = RqLq , 



(2) 



for given values of Rq and Lq [hereafter we will denote 
the pair of parameters {Ro,Lo) by O; we further assume 
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that Rq <C Lq]. For the parabohc shape assumed in this 
paper, the smallest cross-sectional radius is determined 
for any given value of Lq < L < BLq from Eqs. ([l]) and 
(0) as 



Rq 

ao = — 



_l + (3o:^_5) 



1/2 



(3) 



i.e., tto = Rq for L ~ Lq, and ao = (i.e., breakage of 
the wire) for L = 5Lq. 



2. Shell Correction Method 
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The Shell Correction Method we employ is 
the LDA theory. In the Shell Correction MethodE 
(SCM), the total LDA energy, Et{L,0), for any config- 
uration of the wire (specified by L and O) is separated 



Et{L,0) = E{L,0) + AE,i,{L,0) , 



(4) 



where E{L,0) varies smoothly as a function of the sys- 
tem size (L) while AEshiL,0) varies in an oscillatory 
manner with i, as a result of the quantization of the 
electronic states. Ai?sh(riw^) is usually jjpJ led a shell 
correction in the nuclearn3c3 and clusterEHlil literature; 
we continue to use here the same terminology with the 
understanding that the electronic levels in the nanowire 
form subbands, which are the analog of electronic shells 
in clusters where the size of the system is usually given 
by specifying the number of atoms TV. The SCM method, 
which has been shown t©|A4eld results in excellent agree- 
ment with |£pqperimenta3l23'E3 and self-consistent LDA 
calculationsSO for a number of cluster systems, is equiv- 
alent to a Harris functional (i?harris) approximation to 
the Kohn-Sham LDA with the input density, p™, ob- 
tained through variational minimization of an extended 
Thomas- Fermi (ETF) energy functional, Ebtf[p]- 

The Harris functional is given by the following expres- 



-Eharris [p'"] — -E/ + ^ t 
1 



1=1 



FH[p-(r)]-t-14c[p'"(r)] j>p'"(r)dr 
£:.e[p'"(r)]dr , 



(5) 



where Vh is the Hartree (electronic) repulsive poten- 
tial, Ej is the repulsive electrostatic energy of the ions, 
and Ey^c p, / £xc[p]dr is the exchange-correlation (xc) 
functional^ [the corresponding xc potential is given as 
V^c(r) = SExc[p]/Sp{r)]. e°"* are the eigenvalues (non- 
self-consistent) of the single-particle Hamiltonian, 



with the mean-field potential given by 

V,n[p"\r)] = Vh[p"\v)] + V.,[p^^iv)] + Viir) , (7) 

Vi (r) being the attractive potential between the electrons 
and ions. 

In electronic structure calculations where the corpus- 
cular nature of the ions is included (i.e., all-electron or 
pseudo-potential calculations), p™ may be taken as a su- 
perposition of atomic-site densities. In the case of jellium 
calculations, we have showio'Ell that an accurate approx- 
imation to the KS-LDA total energy is obtained by using 
the Harris functional with the input density, p'", in Eq. 
(^ evaluated from a variational Extended-Thomas- Fermi 
(ETF)-LDA calculation. 

The ETF-LDA energy functional, Ebtf [p] , is obtained 
by replacing the kinetic energy term, T[p\, in the usual 
LDA functional, namely in the expression, 

Ei^uaIp] ^ T[p] 

+ J {^lvH[pir)]+Viir)^p{r)dr 



£^M:^)]dr + Ei 



(8) 



by the ETF kiBSitic energy, given to the 4th-order gradi- 
ents as followsJ23 

2me 2rne f 

-^7etf[pJ = / tETFlPldr = 

^2 



3 ^ p 



dr 



/3 



(9) 



The optimal ETF-LDA total energy is then obtained 
by minimization of i^ETrlp] with respect to the den- 
sity. In our calculations, we use for the trial densities 
parametrized profiles p(r; {"fi}) with {7^} as variational 
parameters (the ETF-LDA optimal density is denoted as 
p). The single-particle eigenvalues, {e°"*}, in Eq. (||) 
are obtained then as the solutions to the single-particle 
Hamiltonian of Eq. (||) with Vin replaced by Vetf [given 
by Eq. (0) with /9'"(r) replaced by p(r)]. Hereafter, these 
single-particle eigenvalues will be denoted by {ci}. 

In our approach, the smooth contribution in the sepa- 
ration (||) of the total energy is given by -Eetf[p], while 
the shell correction, A_Esh, is simply the difference 

AEsh= i^harris[p] - EeTf[p\ 
occ p 

= J2^^- I p{r)VETF{r)dr ~ Tetf[p\ ■ (10) 



H = — 



2me 



(6) 
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3. Adiabatic Assumption 

The volume density of the positive background is given 
by p+ = 3/(47rrj), where is the Wigner-Seitz radius 
characteristic to the material, and thus the number of 
positive charges in the constriction is 



(11) 



Since the nanowire contains a constricted region of 
variable cross-sectional radius a{z) [see Eq. (||)], we de- 
fine a linear (i.e., density per unit length of the nanowire) 
background density p^{z;L,0) = 3a^(z)/(4r^), which 
when integrated over the total length of the wire yields 
N'^{0) [see Eq. (|lTl)]. Correspondingly, the variational 
electronic volume density p{x:L,0) = p(r, z; L,0), and 
in our calculations it takes the form, 



p(r,z;L,0) 



Po{z) 



1 + cxp 



/ r-ro{z) \ 



(12) 



with Pq{z), a{z), and j{z) as z-dependent variational pa- 
rameters. In the ETF calculation, p is determined vari- 
ationally at a given z as the one associated with a uni- 
form cylinder of radius a{z) (adiabatic assumption), un- 
der the normalization condition for local charge neutral- 
ity, namely, 2tt J dr[rp{r, z; L,0)] =pt{z;L,0) [which 
fixes the 4th parameter ro{z) in Eq. (|l^)]. The optimized 
p allows then calculation of the smooth contribution for 
any length of the constriction, E{L, O) = E-etf{L, O) in 
Eq. (I). 

The calculation of the shell-correction term, 
AEsh{L,0), in Eq. (^ proceeds by evaluating first the 
density of states in the nanowire. Assuming an adiabatic 
separation of the |iifaatll| transverse and the "slow" longi- 
tudinal variables jljHa the electronic wave functions in 
the classically allowed regions may be written as 



,(r,(/.,z;L,0)a 7/;„™(r;z,L,0)e""* 



xe 



i J' dz'k'Y^(z'\e,L,0) 



(13) 



where fc"™ is the local wave number along the axial (z) 
direction of the nanowire 



fc^'"(z,e;L,0) 



2mp 



1/2 



(14) 



and 'tnm is the (transverse) local eigenvalue spectrum at 
z. To calculate this spectrum for a wire of a configuration 
specified by (L, O) for any value of z, the eigenvalues of 
a cylindrical wire with a (uniform) radius a(z) are calcu- 
lated from the two-dimensional Schrodinger equation 



2mp 



1 d 

r dr 



m 



^ + V^T^{r-z,L,0)ij 



= e 



'„m(z,L,0)V' 



(15) 



The linear (per unit length), one-dimensional density 
of states at z, Di{z,e;L,0), is given by 

2 ^ dkl"'iz,e;L,0) 



Di{z,e;L,0)=-y2- 

7T ^ — ' 

nm 

xe[e-?„™(z;L,0)] , (16) 



TT ^ — ' de 

nm 



where spin degeneracy has been included, and O is the 
Heaviside step function. 

(14), we obtain 



From Eq. 



Di{z,e;L,0) 



2r71e 

1^ 



1/2 



5][e-?„,^(z;L,0)]-^/2 

nm 

xe[e-?„„,(z;L,0)] . (17) 



We may now define an integrated density of states in 
the constriction 



rL/2 

D{e;L,0)^ / dzDi{z,e; L,0) 

J-L/2 



(18) 



The total number of states up to energy e in the con- 
stricted region of the wire is given by 



N-{e;L,0) = f de'D{e'-L,0) 
Jo 



2 f^/^ [2^^ Z 



J-L/2 V h 

xe[e - e„„,(z;L,e')] . 



(19) 



Since the total number of electrons in the constricted 
region is N^iO) [see Eq. (pT|)], the Fermi energy, 
epiLjO), for a wire with a configuration specified by 
(L, O) is given from Eq. (|l9|), i.e.. 



N-{eF;L,0) = N+{0) 



(20) 

Using the above and Eq. (|o|), the shell-correction 
term, 

AS,h(L, O) = Gharris [p] L, O] - E^TF [p] L, O] , (21) 

may be calculated as 

eF(L,0) 



AEsh{L,0) 

cL/2 



de[eD{e;L,0)] 



dz / drrp{r,z;L,0)VETF{r,z;L,0) 

L/2 Jo 

/L/2 /"Oo 
dz / drrtETF [p(r, z; L, O)] , (22) 
-L/2 Jo 

where Vetf is the ETF potential (Hartree, exchange- 
correlation, and electron attraction to the positive back- 
ground) and t-ETF is the volume density of the ETF 
kinetic-energy functional [see Eq. (^] . 
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In actual calculations, we invert the order of integra- 
tion in the first term of Eq. (p2|), which then takes the 
form 

— / V[ej. + 2e„„(z;L,0)] 

^\l^[eF-enm{z;L,0)]Q[ep~e^^{z-L,0)] . (23) 

Note that Eq. (|2^) implies a common Fermi level for 
the whole constriction for a given L (i.e., ep is not a 
local property). Therefore, Eq. ( p3| ) is not equivalent 
to integration of the corresponding uniform wire result 
derived by us in Rcf. |l^ over the z-coordinate, since there 
ep varies with the wire's cross-sectional radius. 

Having calculated the smooth and shell-correction con- 
tributions to the total energy, as a function of L, the to- 
tal elongation force may be evaluated as the derivative of 
the total energy with respect to L, i.e., Ft = —dET/dL, 
and the contributions to it from the ^mooth and shell- 
correction terms are given by F = —dE/dL, and AFsh — 
-dAE^i,/dL. 

B. Results 

In this section, we report results for the elongation of 
a Sodium nanowire (rg = 4 a.u.), starting with an initial 
cylindrical constriction of length Lq — 80 a.u. and radius 
Rq = 25 a.u. 

In Fig. 2 , we show electronic-potential profiles, 
VETFif'TO-iz), L,0], for a particular constriction with 
AL/Lq ^ 1.125 (AL = L - Lq). We display here the 
potential profiles calculated at the narrowmost part of 
the constriction [qq = a(0) — 12.62 a.u.] and at its end 
[i.e., at a{L/2) = Rq] (in this paper, all the subsequent 
numerical results we will discuss relate to constrictions 
with the same set of O-parameters, namely, Lq = 80 a.u. 
and Rq — 25 a.u.). We found that for other values of 
z (i.e., for < [zl < L/2), the potential assumes pro- 
files intermediate between the two profiles shown here, 
namely, the depth of the potential well remains practi- 
cally unaltered, while its width follows the enlargement 
of the jellium-background radius a{z), from qq to a{L/2). 
From the three components, Vh, Vj, and V^c, which con- 
tribute to the total Vetf [see Eq. (^], we found that for 
all these potential profiles, calculated for different val- 
ues of z along the constriction, the xc contribution is the 
dominant one, amounting to approx. —5.4 eV, while the 
total electrostatic contribution, Vh + Vj, is much smaller, 
resulting in a characteristic "winebottleU, profile familiar 
from LDA studies of spherical clusters.EEl 

Fig. 2 also displays the transverse local eigenvalues 
Enm associated with the two potential profiles. Natu- 
rally, a wider potential profile yields a larger number of 
such eigenvalues below the Fermi level. 

To illustrate the nature of the electronic spectrum in 
the nanowires, and its dependencies on the characteristics 



of the wire, i.e., shape and length, we show in Fig. 3(a) 
densities of states, D{e; L,0), calculated for a variable- 
shaped wire for two wire lengths (and consequently two 
minimal cross-sectional radii of the constricted region). 
The density of states for a uniform wire with a radius 
equal to that of the unconstricted region of the variable- 
shaped ones [shown in Fig. 3(a)] is displayed in Fig. 3(b). 

Two "classes" of features are noted for the variable- 
shaped wires: (i) those associated with the narrowmost 
constricted region (marked by numbers) whose radius, 
uq, varies upon elongation, and (ii) those associated with 
the maximal radius of the constriction (and with the 
unconstricted part of the wire) which remains constant 
throughout the elongation of the wire. Identification of 
the latter class of features (several of which are marked by 
arrows) is facilitated through comparison with the den- 
sity of states for the corresponding uniform wire [Fig. 
3(b)]. We observe here that, for the broader (and thus 
shorter) wire [lower curve in Fig. 3(a)], six of the fea- 
tures (peaks) in the density of states coming from the 
spectrum of transverse energy levels at the narrowmost 
region of the constriction are located below the Fermi 
level, ep [all the peaks in the density of states occur at 
the energies of the transverse levels; e.g., compare the 
location of the peaks in the lower curve in Fig. 3(a) with 
the corresponding spectrum on the left side of Fig. 2] . 

On the other hand, for the much narrower (and thus 
longer) constricted wire, only one of these peaks is below 
ep [see upper curve in Fig. 3(a)]. When plotting the den- 
sity of states at ep versus the elongation (or equivalently 
the minimal radius of the constricted region) , these vari- 
ations lead to an oscillatory pattern, as peaks in the den- 
sity of states are shifted above the Fermi level, one after 
the other as the wire is being elongated. These variations 
are also portrayed in the energetics of the wire (shown 
in Fig. 4), and in the stepwise behavior of the quantized 
conductance through the wire versus length (see Fig. 5 
below) . 

From Fig. 4, we observe that the magnitude of the 
smooth ETF contribution, E, to the total energy, Et, of 
the wire is dominant, with the shell-correction contribu- 
tion, AEsh, exhibiting an oscillatory pattern, with local 
minima at a set of wire lengths (and correspondingly a set 
of minimal cross-sectional radii) which we term "magic 
wire configurations" (MWC's), i.e., wire configurations 
with enhanced energetic stability. When added to the 
smooth contribution, these shell-correction features lead 
to local minima of the total energy toward the end of the 
elongation (and consequently, narrowing) process, while 
for thicker wires (i.e., AL/Lq < 2.5 in Fig. 4) they are 
expressed as infiection points of the total energy (in this 
context, see the total-force curve. Ft, in Fig. 5, where 
the local minimum in Ft corresponds to the point with 
Ft = marked by an arrow). 

We note here that the occurrence of local minima in 
the total energy results from a balance between AEgh 
and E, with the latter increasing (that is acquiring less 



5 



negative values) as the constriction elongates due to the 
increasing contribution from the surface of the constric- 
tion. Comparison of the magnitudes of the shell cor- 
rections in a variable-shaped wire and in a uniform one 
[i.e., one with /(z) = 1 in Eq. (|^), whose case was dis- 
cussed in Ref. Q shows that the amplitudes of the oscil- 
lations in the latter case are much larger (over an order of 
magnitude). The reason for this difference is that in the 
constant-radius wire the quantization into the transverse 
subbands is uniform along the wire, while in the variable- 
shaped case the subband spectrum is different in various 
parts of the constriction. While the oscillatory pattern 
is dominated by the spectrum at the narrowmost region 
(see also Section III below), the amplitudes are influ- 
enced by the transverse-mode spectra from other parts 
of the constriction. Consequently, the number of local 
minima in the total energy, Et, (and thus the number 
of wire configurations, i.e., lengths, for which the total 
force. Ft, vanishes) is larger for a uniform wire than for 
a variable-shaped one. Additionally, we suggest that for 
materials with relatively smaller surface energies a larger 
number of local minima may occur. 

From the total energy, and the smooth and shell- 
correction contributions to it, we obtain the total "elon- 
gation force" (EF), Ft, and the corresponding compo- 
nents of it, F and A_Fsh- These results are displayed in 
Fig. 5, along with the conductance of the wire evaluated, 
in the adiabatic approximation (i.e., no mode mixingtS) 
and neglecting tunneling effects (assuming unit transmis- 
sion coefficients for alLthe conducting modes), using the 
Landauer expression ,E3E3 

G(L, 0)^goY, ^i^F - enmiz = 0; L, O)] , (24) 

nm 

where go = 2e^ /h, and the spectrum of the transverse 
modes is evaluated (for each constriction length) at the 
narrowmost part of the constriction, z = 0. Tunneling 
contributions (see e.g., Ref. mode-mixing and non- 
adiabaticity may affect the sharpness of the conductance 
steps, and/or introduce some interference related fea- 
tures, particularly near the transitions between the con- 
ductance plateaus. These effects, which can be indjuded 
in more elaborate evaluations of the conductance,&Ell 
do not modify the conclusions of our study. 

Also included in this figure is a plot describing the 
variation of the minimal cross-sectional radius aq with 
the length of the constriction [see Eq. (^)] . 

As evident from Fig. 5, the oscillations in the force re- 
sulting from the shell-correction contributions are promi- 
nent. In AFsh, we observe that the locations of the ze- 
roes of the force situated at the right of the force max- 
ima occur for values of AL/Lq which coincide with the 
locations of local minima in the shell-correction contri- 
bution to the energy of the wire (i.e, for a sequence of 
minimal cross-sectional radii corresponding to MWC's). 
In the total force. Ft, only one of these points (where 
Ft — 0) remains [i.e., the one corresponding to the lo- 



cal minimum in the total energy towards the end of the 
elongation process (see Fig. 4)], for the reasons discussed 
above in connection with the energetics of the wire. Nev- 
ertheless, the oscillations in the total force correlate well 
with those in the total energy of the wire, which as dis- 
cussed above originate from the subband spectrum at 
the narrowmost part of the constriction (see also section 
III). Also, the locations of the local maxima in the total 
force correlate with the stepwise variations in the conduc- 
tance signifying the sequential decrease in the number of 
transverse subbands (calculated at the narrowmost sec- 
tion of the wire) below tp (i-e., conducting channels) as 
the constricted part of the wire elongates (and thus nar- 
rows). Additionally, we note that the magnitude of the 
total force is comparable to measured ones (i.e., in the 
nanonewton range). The magnitude of the total force in 
sodium nanowires (not measured to date) is expe^fted to 
be smaller than that found for gold nanowiresJ^O due 
mainly to differences in the electron densities and surface 
energies of the materials. 

III. SEMICLASSICAL ANALYSIS 

As discussed above, the total energy of the wire is 
characterized by local minima and inflection points oc- 
curing for a set of wire lengths, or equivalently a set 
of minimal cross-sectional radii of the constriction, and 
are reflected in the oscillatory patterns of the elongation 
force. These features correspond to the oscillatory shell- 
correction contributions and originate from the spectrum 
of transverse modes at the narrowmost part of the con- 
striction. Moreover, these patterns correlate with the 
locations of the quantized conductance steps, which are 
determined by the transverse-mode spectrum at the nar- 
rowmost region (i.e., the number of conducting modes 
below €f, and their degeneracies). 

To further investigate the origins of these correlations, 
we present in this section a semiclassical analysis of the 
density of states, energetics, forces, and conductance in a 
free-electron nanowire modeled via an infinite confining 
potential on the surface of the wire. As in the above 
(see Fig. 1), we model the constricted region of the wire 
as a section with a slowly (adiabatically) varying shape. 
Dividing the constriction into thin cylindrical slices, the 
solution of the Schrodinger equation for each slice is of 
the form, 

V; = ytJ„(«;r)e™'^e*f^^ , (25) 

where ^ is a normalization constant, pj_ is the electron 
momentum along the axis of the wire, Jm^Kr) is the 
Bessel function of order m, and n = (2mee — p\y^'^ /h. 

Consider first a uniform cylindrical wire with a con- 
stant cross-sectional radius a. With the infinite wall 
boundary condition assumed here, the single-particle 
electronic energy levels in the wire are expressed in terms 
of the roots of the Bessel functions, 7„m, as 
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2 

nm 



+ 



Pi 



2mea? 2m 



(26) 



Here we remark that in the semiclassical approxima- 
tion the electron performs a compHcated trajectory inside 
the wire. AU the semiclassical trajectories are tangent to 
the caustie surfaces of a set of concentric cylinders inside 
the wircca Quantization of the electronic states leads to 
selection of only a certain subset of trajectories associated 
with a ceratin set of radii of the caustic surfaces, cor- 
responding to allowed values of the azimuthal quantum 
numbers, m, i.e, Kr,„ — m; this description is closely re- 
lated to the semiclassical periodic orbit theory.H In the 
course of-developing semiclassical methods, Keller and 
RubinoweSl havje-demonstrated that the Debye asymp- 
totic expansioneil of the Bessel functions (1 ^ to < Kr) 
provides an accurate approximation to the eigenfunction 
Jminr), i.e., 

^™W-(-)'/'(«V-to2)-V4 

vr 

(K^r^ - TO^)^/^ - TOarccos ( — ) 4- - . (27) 

\Kr/ 4 J 



xsm 



This approximation is valid in the region between the 
caustic cylindrical surface and the boundary surface of 
the wire; in the region inside the caustic surface (to > kt) 
the solution decays exponentially. In this approxima- 
tion, the equation for the asymptotic values of the Bessel- 
fimction zeroes has the form, 



(7. 



2 

nm 



,2n1/2 



TO arccos 



In 



= TT n — 



(28) 



First we calculate the density of states whose evalu- 
ation involves, after integration over p±, double sums 
over the quantum numbers n and to; n = 1, 2, 
TO = 0, ±1, ±2, ... [see Eq. (|l7|)]. Applying sequentially 
the Poisson summation formula to both sums and sep- 
arating the oscillatory terms (note that in our semiclas- 
sical approximation Ka 3> 1) in complete analogy with 
Refs. ||,|6[ we obtain for the density of states (per unit 
length) , 



nrie) 



iraea 



oo M/2 ^ 

X > > — si 

M=2 Q=l 



sml I cos 



2MKa sinl 



ttQ 
M 



■kM' 



3/4 1/4 



— T — 

1/4 .4^ Ml 



M=l 



/2 



2TTMKa + - 



(29) 



where Ca = /{2mea'^), and K is the electron wave vec- 
tor. The two terms in Eq. ( p9| ) correspond to the contri- 
bution from the point where the phase is stationary and 
from the end-points in the sum (integral) over m (see 
discussion in Ref. |4^). While the second oscillatory term 
in Eq. (E9h has a smaller amplitude than the first one 



[by a factor of {Kay^'^], it corresponds to an important 
class of electronic states, with to « Ka, localized near 
the suxface of the wire (the so-called whispering gallery 
statescj). 

Until now we discussed a uniform wire with a constant 
cross-sectional radius. In a wire with a variable shape, 
the cross-sectional radii depend on z, as discussed in con- 
nection with Eq . (|l| ). Substituting the z-dependence of 
the radii in Eq. (^9D, i.e., replacing a by a(z), we need to 
perform an integration over z [see Eq. (p^j. This inte- 
gration involves evaluation of integrals of the form, 



L/2 

I = 3t I .g(z)e'"-^'^Wdz , 

-L/2 



(30) 



where for the first term in Eq. ( p9| ) g{z) — a{z) and a — 
2Msm{TTQ/M), and for the second one g{z) = \/a{z) 
and a = 2-kM . The fast oscillatory character of the ex- 
ponential factor [i.e., Ka{z) ^ 1 for all z] relative to 
the slow variation of ^(zV-allows us to use the standard 
stationary phase method,cZI obtaining 



27r 



1/2 



aKa"{0) 
2 

' aKa'{L/2) 



g(0) 3?{exp[iai4:a(0) + i7r/4]} 
g(L/2) 5R{-iexp[mifa(L/2)]} , (31) 



where z = is the stationary (extremum) point, the sec- 
ond term is the contribution from the end-points of the 
integral, and primes denote differentiation with respect 
to z. Using the above, and after simple algebraic manip- 
ulations, we obtain for the oscillatory part of the density 
of states, 



D°^'=(e) 
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sm 
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1 2mML 2) 

--rT^—!i sm 

M2 n^Ka'{L/2) 
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M 
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1 
' 2 
ttM 
~2~ 
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1/2 ^ 


1m 


Ka"{0) 


a(0) 



3/4 



cos[27rMi^a(0)] 



tM^/^ Ka'{L/2)a{L/2) \ 



f 2mea^{L/2) \ 



3/4 



X COS 



2-KMKa{L/2) + - 



}■ 



(32) 



The density of states of the wire contains oscillatory con- 
tributions from the narrowmost cross-section of the wire 
[first and third term in Eq. (|32|)] and from the wire's-end 
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cross-sections (second and fourth terms) . The amphtudes 
of the latter oscillations is smaller. 

Having obtained the expression for the oscillatory part 
of the density of states, we can calculate now the semi- 
classical approximation to the grand-canonical thermo- 
dynamic potential O [see Appendix A; at zero temper- 
ature, n = /(e — eF)D{e)de]. Restricting ourselves for 
brevity to the largest contribution [that is, to the first 
term in Eq. ( ^ ) corresponding to the main contribution 
from the narrowmost part of the wire], we get for the 
oscillatory part of fi, 



2eF 



oo M/2 ^ 

V^kFa"iO) a{0) £ MV^ 



sm 



-3/2 



(^g/Af) 



2MkFa{0) sin(7rg/M) + f ( + ^ 



(33) 



From this expression, the oscillating part of the force 
as a function of the length of the constricted region [i.e., 
a(0) in general depends on L, see e.g., Eq. (||)] is given 
by 



dn""" da{0) 



daiO) dL ' 
which upon substitution of (|3^) yields, 

4eFk]/^[da{0)/dL] ^ ^ 1 

M=2 Q=l 



(34) 



X cos 



y/na"{0) a(0) 

2MkFa{0) sin(7rg/Af) 



-i/2(7rg/Af) 



M - 



(35) 



The expression for the conductance of the wire fol- 
lowing the Landauer formula involves evaluation of the 
number of transverse states in the narrowmost part of 
the wire. Following Ref. n§, 



G: 



2e^ 
h 

oc. M/2 



[kFa{Q)]' 



1 - 



1 



kpaiO) [fcFa(0)]3/2 



M=2 Q=l 



3ini/2(7rg/M) 



2MkFa{0) sininQ/M) 



M 



(36) 



which can be expressed as a function of the length of 
the constriction [see e.g., Eq. (||); we remark here that 
our semiclassical treatment is valid for any adiabatic wire 
shape] . 

The non-oscillating contribution (coming from the fir|St 
two terms in the curly brackets) describes the-Sharvinua 
conductance of the constriction and the WeyO semiclas- 
sical corrections, and the third term describes conduc- 
tance quantum oscillations as a function of a(0). From a 



comparison of the expression for the oscillatory contribu- 
tion to the force [Eq. (|35|)] with the oscillatory contribu- 
tion to the conductance [Eq. (^)], the direct correlation 
between the two is immediately evident, and both de- 
pend on the spectrum of transverse modes (conducting 
channels) at the narrowmost part of the wire. This is 
in agreement with the results shown in Fig. 5 obtained 
through the LDA-SCM method. 



IV. CONCLUSIONS AND DISCUSSION 

In this paper, we extended our investigation^ of en- 
ergetics, conductance, and mesoscopic forces in a jel- 
lium modelled nanowire (sodium) using the local-density- 
functional-based shell correction method to variable- 
shaped wires, i.e., containing a constricted region mod- 
eled here by a parabolic dependence of the cross-sectional 
radii in the constriction on z (see Fig. 1). The results 
shown above, particularly, the oscillations in the total 
energy of the wire as a function of the length of the 
variable-shaped constricted region (and correspondingly 
its narrowmost width), the consequent oscillations in 
the elongation force, the corresponding discrete sequence 
of magic wire configurations, and the direct correlation 
between these oscillations and the stepwise quantized 
conductance of the nanowires, originate from quantiza- 
tion of the electronic states (i.e., formation of subbands) 
due to the reduced lateral (transverse) dimension of the 
nanowires. These results are in correspondence with our 
earlier LDA-SCM. investigation of jellium-modeled uni- 
form nanowires. E2I Moreover, in the current study of a 
wire with a variable (adiabatic) shaped constriction, we 
found that the oscillatory behavior of the energetic and 
transport properties are governed by the subband quan- 
tization spectrum (termed here electronic shells) at the 
narrowmost part of the constriction. This characteristic 
is supported and corroborated by our semiclassical anal- 
ysis (section III). 

We reiterate here that such oscillatory behavior, as 
well as the appearance of "magic numbers" and "magic 
configurations" of enhanced stability, are a general char- 
acteristic of finite-size fermionic systems and are in direct 
analogy with thospiound in simple- metai clusters (as well 
as in 3 He clusteraHa and atomic nucl^ 
tronic shell effects on the, energetics 
recently shape dynamicsEil of jellium modelled clusters 
driven by forces obtained from shell-corrected energet- 
ics) have been studied for over a decade. 

While these calculations provide a useful and instruc- 
tive framework, we remark that they are not a substi- 
tute for theories where the atomistic iiatiirg|-^Lnd spe- 
cific atomic arrangements are includedliJiij'ErEl in eval- 
uation of the energetics (and dynamics) of these sys- 
tems (see in particular Refs. pUq, where first-principles 




^ where elec- 
£3 (and most 



molecular-dynamics simulations of electronic spectra, ge- 
ometrical structure, atomic dynamics, electronic trans- 
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port and fluctuations in sodium nanowires have been dis- 
cussed). 

Indeed, |— the, atomistic structural characteristics of 
nanowirealjiij (including the occurrence jo£ cluster- 
derived structures of particular geometriesBB) , which 
may be observed through the use of high resolution 
microscopy^a influence the electronic spectrum and 
transport characteristics, as well as the energetics of 
nanowires and their mechanical properties and response 
mechanisms. In particular, the mechanical response of 
materials involves structural changes through displace- 
ment and discrete rearrangenent of the atoms. The mech- 
anisms, pathways, and rates of such structural transfor- 
mations are dependent on the arrangements and coordi- 
nations of atoms, the magnitude of structural transfor- 
mation barriers, and the local shape of the wire, as well as 
possible dependency on the history of the material and 
the conditions of the experiment (i.e., fast versus slow 
extensions). Further evidence for the discrete atomistic 
nature of the structural transformations is provided by 
the shape of the force variations (compare the calculated 
Fig. 3(b) in Ref. |l| and Fig. 3 in Ref. |l| with the mea- 
surements shown in Figs. 1 and 2 in Ref. ^), and the 
interlayer spacing period of the force oscillations when 
the wire narrows. While such issues are not addressed 
by models which do not include the atomistic nature of 
the material, the mcsoscopic (in a sense universal) phe- 
nomena described by our model are of interest, and may 
guide further research in the area of finite-size systems in 
the nanoscale regime. Such further investigations include 
the occurrence of magic configurations (i.e., sequences of 
enhanced stability specified by number of particles, size, 
thickness or shape) in clusters, dots, wires, and thin films 
of normal, as well as superconducting metals, and the ef- 
fect of magnetic fields which can influence the energetics 
in such systems (e.g., leading to magnetostriction effects) 
through variations of the subband spectra, in jajjialogy 
with magnetotransport phenomena in nanowirea23. 

Several directions for improving the model (while re- 
maining within a jellium framework) are possible. These 
include: (i) consideration of more complex shapes. For 
example, in our current model the elongation is dis- 
tributed over the entire constriction throughout the pro- 
cess, while a more realistic description should include a 
gradual concentration of the elongation, and consequent 
shape variation, to the narrower part of the constriatioB, 
as found through molecular-dynamics simulations ?13^Ell 
(ii) use of a stabilized-jellium descriptioij£3 of the ener- 
getics of the nanowire in order to give it certain elements 
of mechanical stability. In this context, note also that 
from the total energy shown in Fig. 4(c), and the corre- 
sponding total force [Fig. 5(c)], it is evident that in our 
current model, except for the region of large elongation 
close to the breaking point (i.e., AL/Lq > 2.5), the wire 
is unstable against spontaneous collapse (that is short- 
ening), i.e., there are no energetic barriers against such, 
process, while both experimentsEj and MD simulationsE^ 
show that compression of such wires requires the applica- 



tion of an external force. Improvements of the model in 
these directions are moEjt|desirable in light of tiie afore- 
mentioned experimenta]E3 and MD-simulationsll2l obser- 
vations that the total oscillating forces for elongation and 
compression of nanowires are of opposite signs (i.e., neg- 
ative and positive, respectively), while our current (equi- 
librium model) is limited to certain aspects of the tensile 
part of an elongation-compression cycle; (iii) inclusion of 
bias voltage effects in calculations of the energetics and 
conductance of nanowires .oEj While such effects may be 
expected to have little influence (particularly on the en- 
ergetics) at small voltages, they could be of signiflcance 
at larger ones. Work in these directions is in progress in 
our laboratory. 
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APPENDIX: A 

In this Appendix, we discuss briefly a semiclassical 
treatment of temperature effects on the oscillatory be- 
havior of the force and conductance in nanowires. The 
grand-canonical thermodynamic potential at finite tem- 
perature, T, is given by 



-k, 



.T^ln 



1 + exp 



(Al) 



where i denotes (n, m,p±), and fi is the chemical poten- 
tial. 



From Eq. (Al), the finite temperature expressions for 
fiosc^ posc^ g^j^^ Qosc (jjffgj. from those given for the zero- 
temperature limit in Eqs. (^, (p5|), and (36), respec- 
tively, by a multiplicative factorjjn the sums of these 
equations. This factor is given by^^ 



*(^mq) 



X 



MQ 



sinh( Amq) 



where 



Xmq — 



2TrMkBTa{0) sitl{ttQ/M) 



(A2a) 



(A2b) 



with vp being the Fermi velocity. For T — 0, 'I'(.x) — 1. 

N ote that the temperature dependence given in Eq. 
( |A2| ) is valid for systems with kpaiO) ^ 1, and leads to 
reduction of the oscillation amplitudes when 2TrAIkBT > 
Ae, where Ae = hvp /[a{0) sm{TTQ/M)] is an effective 
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energy-level spacing of the electrons contributing to the 
oscillatory parts of the thermodynamic potential, force, 
and conductance. 
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FIG. 3. Densities of states for: (a) two configurations 
of the variable-shaped wire, one (lower curve) with elonga- 
tion AL/Lo — 1.125 and narrowmost radius ao = 12.62 
a.u. (potential profiles and local transverse spectra for this 
case are displayed in Fig. 2), and the other with elongation 
AL/Lo = 2.75 and narrowmost radius ao = 4.57 a.u. (up- 
per curve whose y-axis is shown on the right), (b) a uni- 
form- in-shape wire with AL/Lo = and ao = -Ro = 25 
a.u. For all cases, Lq — 80 a.u. and Tio = 25 a.u. The 
vertical dashed lines denote the corresponding Fermi lev- 
els. The Fermi level of the constrictions, which for the uni- 
form-in-shape wire is —2.82 eV, varies only by 0.05 eV for 
all the elongations down to the breakage point. In (a) the 
numbered peaks correspond to the locations of the transverse 
energy levels in the narrowmost part of the constriction [e.g., 
compare the lower curve in (a) with the spectrum shown on 
the left axis of Fig. 2]. The arrows indicate the locations 
of some of the transverse energy levels at the end-points of 
the constriction, coinciding with corresponding peaks in the 
spectrum of the uniform-in-shape wire shown in (b) . 



FIG. 4. Energies (in eV units) of a variable-shaped sodium 
nanowire, plotted versus the relative elongation AL/Lq. The 
initial parameters O are Lo = 80 a.u. and Ro = 25 a.u. The 
smooth, ETF contribution (E), the shell correction (AEsh), 
and the total energy (Et) are displayed in (a), (b), and (c), 
respectively. 



FIG. 1. Schematic drawing of the jellium-background 
of a variable-shaped nanowire. The (cylindrical) symme- 
try axis of the wire is along the z-axis, with a con- 
stricted region (— L/2 < z < L/2) described by a depen- 
dence of the cross-sectional radii a(z) on z [see Eq. (^)]. 
Ro — a{z — ±L/2) is the radius in the uniform part of the 
wire outside the constriction. 

FIG. 2. Potential profiles at the narrowmost {z = 0; left 
curve) and end-points {z — ±L/2; right curve) of a constric- 
tion with an elongation AL/Lq = 1.125 (whose O-parameters 
are Lo = 80 a.u. and Ro = 25 a.u.). The local trans- 
verse-mode spectra, e„m, associated with these two profiles 
are also displayed along the left and right y-axes. The dashed 
line indicates the Fermi level. For the {z = 0)-profile, the 
spectrum is labeled with the corresponding n, m local trans- 
verse eigenvalues (n denotes the number of radial nodes plus 
one, and m denotes the azimuthal angular-momentum quan- 
tum number). The same spectrum is numbered sequentially 
(in parentheses) to facilitate comparison between the location 
of the energy levels and the numbered peaks in the density of 
states given in the lower curve in Fig. 3(a). Energies in units 
of eV, and lengths in a.u. 



FIG. 5. (a-c): The smooth, ETF contribution to the force 
(F), shell-correction force (AFsh), and total force (Lr), corre- 
sponding to the energies shown in Fig. 4(a-c), plotted versus 
the relative elongation AL/Lo. The arrow in (c) indicates 
the point Ft ~ corresponding to the local minimum in the 
total energy shown in Fig. 4(c). The dashed lines indicate 
the zeroes of the y-axes. Forces in units of nanonewtons. (d): 
The conductance, G, for the variable-shaped wire in units of 
go = 2e^Jh, plotted versus AL/Lq, evaluated as described 
by Eq. (PJ). (e): The variation of the cross-sectional radius 
(in a.u.) of the narrowmost part of the constriction, plotted 
versus AL/Lq [see Eq. (^), with Lq = 80 a.u. and Ro = 25 
a.u.] 
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